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1 Introduction 

In this paper we revisit some aspects of a rather classical analytical model occurring 
in many physical settings, the following stationary nonlinear problem: 



§F+iVH = 0, 






in Q, 


u = 0, 






on dfl 


with N(u) := -Lu - 


- Xu - 


-sF(u), 





(1) 



where fi is an open domain in R n with piecewise smooth boundary, u e H := 
H (J?; R), -Lu is elliptic and F(u) is a nonlinear operator, with e > a perturba- 
tion parameter. Here -Xu represents a sort of linear soft device, possibly breaking, 
for A > large enough, the elliptic character of -Lu - Xu. Stationary elastic mem- 
branes, reaction-diffusion, Ginzburg-Landau, among others, can be settled into 
this format. 

First at all, we recall that if N'(u) is symmetric with respect to the L 2 scalar 
product, (F'(u) actually) Volterra-Vainberg theorem applies (see [23], U), and 
we can write the energy functional: 

E{u) = £{N{tu),u) L2 dt, (2) 

and obtain a variational formulation for ([T]): dE(u) = 0. Then, we apply to E(u) a 
version of Amann-Zehnder (AZ) exact finite reduction, which is a global Lyapunov- 
Schmidt technique, introduced in Hamiltonian dynamics in [2] and [10] . and in 
fields theory in pQ (see also [5ll6lll8]). ([1]). In such a way, we reduce the original 
variational problem to a finite dimensional system of 'algebraic' equations: the en- 
ergy functional (J2j) is shrunk to a m variables function E : R m ->■ R, whose critical 
points, i.e., dE(x) = 0, correspond exactly to the solutions of (JXJ) . 

We will introduce some appropriate conditions on F in order to prove that the 
reduced functional E(x) is a generating function quadratic at infinity (GFQI), i.e. 
it coincides, or corresponds in a suitable sense (see details in Section [4]), with a 
non-degenerate quadratic form x J Qx outside a compact setjj This is crucial for 
our aim: indeed, proving that the sublevel sets 

Q c = {x : x J Qx < c} , E c = {x : E{x) < c} , 

are diffeomorphic provided \c\ being suitably large, we can apply the 'deformation 
lemma' and claim the existence of a critical point of E corresponding to the critical 
point of Q (0, of course), the so-called rainraax solution, which is in general a 
saddle. This is precisely the case when the linear soft device is strong enough, so 
the original energy functional is not coercive anymore. Actually, E can be very 
different from Q in a compact subset, where many other critical points may arise. 
Nevertheless, the topological complexity of such a compact perturbation E of Q 
cannot wipe out the minmax critical point inherited from the of Q. 

Contemporary literature displays a plenty of existence, multiplicity, bifurcation 
results for nonlinear PDEs like ([1]); the functional setting can be very general, and 



Generating functions quadratic at infinity were originally developed by Chaperon, Siko- 
rav and Viterbo studying Hamiltonian dynamics (see [T|j 8.15 20*21 22 24 ), leading to major 
results in dynamical systems and a new powerful format in symplectic topology. 
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diverse topological technique^ have been applied. Even though the above abstract 
settings allow for very general results of existence and multiplicity, here we propose 
a fairly self contained introduction to AZ and we adopt some stronger hypotheses 
on F, transforming the variational functional ([2]) to a finite GFQI, which finally is 
well arranged towards numerical applications, as e.g. already implemented in [Bj. 
In the last part of the present work (Section [5J, we point out the interplay 
among discrete/continuous and approximate/exact aspects in problem ([1]) , as syn- 
thetically illustrated in the following diagram: 



thermodynamic 
limit N -» oo 



| N - oscillator chain j j oo - elastic string j 

I discrete & approximate I I continuous & exact I 

analogy equivalence | AZ - reduction 



N — oscillator chain 
discrete & 'exact' J 



coordinates 
change 



N — spectral particles 
(abstract) 

discrete & 'exact' 



Fig. 1 In this picture 'exact' means equivalent from the variational point of view of the elastic 
string. 



As a starting point, we recall that, for instance, the elastic string can be mod- 
eled by means of a chain of N beads connected by springs of length h and sending 
N->oo and h -* keeping mass and strength densities constant (see [131 Chapter 
4.]). Here, the continuous system modeled by ([TJ) can be thought as the ther- 
modynamic limit of a large collection of interconnected small particles. Then, by 
applying AZ, we get back to a finite dimensional problem where approximations 
are absent, in the sense that the full nonlinear complexity of the set of equilibria 
is unaltered. However, as described in Section [2] the reduced system is written in 
spectral variables x, therefore most of the mechanical sense of the starting problem 
is encrypted in an abstract representation. 

In order to tailor a physically meaningful interpretation to E(x), we consider 
its leading term, the quadratic part. Surely, it is possible to find a plenty of fi- 
nite linear physical systems whose potential energy corresponds precisely to this 
quadratic part. In Section 16.21 we suggest a suitable linear coordinatqj transfor- 



Conley index[9], Leray-Shauder degree theory, Morse theory I17II12II1T1 . Mountain Pass 
Lemma and Linking theorems [3]. 
3 Actually, we will use the symbol fi for the spectral coordinates instead of x, sec J30K . 
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mation x *+ y = Yx, selecting, among many variationally equivalent systems, the 
one mimicking most faithfully the linear part of the original mechanical problem. 
More precisely, by choosing the above variables y, we pick out the fittest elastic 
chain from the thermodynamic limit sequence converging to the linear part of ([T]) . 
By doing so we provide a further motivation for the AZ reduction coming from 
the microphysics of the matter. This excursion gives in a sense an accomplishment 
to a line of thought started in [19] . 



2 The Amann— Zehnder finite reduction for fields theory 

Let L be an elliptic operator, i.e., 

The scalar product on H, defined as 

«•,•»: HxH^R, 

(u,v) i — >({u,v)):= / -Lu-vdx. 
Jn 

is positive definite, indeed: 



(3) 



{(u,u)):= / (— Lu)udx = - / V • (ua\/u) + / Vu bVm 
Jn Jn Jn 

= - u a\/u-n+ / Vu T aVu > 0. (4) 

Jen ^ - — , — - Jn 
=o =o 

Dirichlet Neumann 

As a result, we will take \\u\\ := {(u,u)) as the norm on H. 
The eigenfunctions associated to the eigenvalues of the elliptic operator L form 
a basis for the space H, 

= A < Ai < A 2 < . . . , 

II A? II = 1- 
For every v e H we can write 

oo oo 

We are able now to define the Green operator of L, g : H -* H, g = (-L) -1 , 

gv = g\ ^VjUj = £ -^-Uj. 

\3-l I J-l A J 

It is clear that 

-Lgv = g(-L)v = v, \/v £ H. 
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For every fixed mtPi, we can consider the following decomposition of H: 

h = r m H © q m H, 

m oo 

V = V m V + <Q m V = ^ VjUj + ^ VjUj . 
j = \ j=m+l 

We will briefly write 

m oo 

7 = 1 j =771+1 

and we will refer to fi e V m H as to the finite head of v, whereas to 77 6 Q m H as to 
the infinite tail of v. 

If we substitute the relation u = gv in the Dirichlet problem ([TJ) : 

- Lu- Xu = eF(u), u\ = (5) 

we obtain 

v- Xgv = eF(gv), 
v=(I-Xgy 1 eF(gv), 

which solutions are the fixed points of the map: 

H — > H, 

v^il-Xg^Figv). 

The map (I - Xg)~ is well defined for A * Xj, j = 0, 1, . . . and can be written as 
(I - Xgy'v = (I-Xgy 1 ({v,}^) = L^A . 

The cut-off decomposition splits the problem into a finite and an infinite part: 

v = ¥ m v + Q m v = fi + 77, (6) 

l i = P m (I-Xg)- 1 sF(g( f i + V )), (7) 

7? = Q m (I- Xg)- 1 eF{g{ l i + T 1 )). (8) 

These equations can also be equivalently rewritten as follows, 

V m N(g(fi + V )) = 0, (finite) (9) 

Q m N(g(^ + V )) = 0. (infinite) (10) 



Under appropriate assumptions on F, (|10p is uniquely solved for 77. Indeed, we can 
prove that 
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Proposition 1 Let F : H -»■ H be Lipschitz, then for every fixed e > and /x £ 
P m i7 there exists m e N suc/i i/iat 



7 7 ^(I-A ff )- 1 £ Q m F(g( M + 7 7 )), 



(11) 



zs a contraction. 



Proof First, we notice that r\ h> g(ii + 77) is a Lipschitz map with constant -j- 

We deduce that (I- A<j) _1 is bounded, indeed 



|(I-A 9 )-\;|| = 



Gfr«) 



< sup 



A, -A 



and clearly sup t— it : = Cj < +00, because the sequence , \ -> 1. By assumption, 
_F is Lipschitz with constant C2, therefore also the map (JTTJ) is Lipschitz: 

|(I - Xg)- 1 eQ m F(g(^ + 771)) - (I - Xg)' 1 eQ) m F (g{ix + 772))! 

pn-n n 

I1 1)1 - TfA 



< 

SCIC2 

A,T1+1 



For sufficiently large m, we conclude that g 5^ 1 2 < 1. n 

Substituting the fixed point 77 = ?7(/i) of map (|11|) into equation (0, we obtain 

/U =(I-A(j)- 1 £P m ^( 9 ( /J + 7]( M ))), M =(/ J l,...,/ Jrn )6R" 1 , (12) 

which is a finite dimensional equivalent formulation for problem ([I]). Indeed, for 
every solution fi of this last equation we can write a solution of ([TJ by means of 
the formula: 

u = g{fi + i){f»)) 
Conversely, to every solution u of ([TJ) , corresponds a solution ft of (|12[) by 

fi = F m (-L)u. 



3 Infinite dimensional and reduced finite dimensional variational 
formulation. 

If we assume the Gateaux derivative N'(u) of the non linear operator considered 
in JT]) to be symmetric with respect to the L -scalar product, i.e. 

(N'(u)h,k) = (N'(u)k,h), Vh,keH, 

(u,v) '■= (u,v) L 2 = / uvdx, 

the Volterra— Vainberg theorem allows us to write an energy functional, 
E:H — -K, 



Tii — >E(u):= / (N(tu),u)dt, 
Jt=o 



(13) 



and the corresponding variational principle dE(u) = is equivalent to the original 
Dirichlet problem ([I]). More precisely, 
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Theorem 1 Every critical point of E is a solution of (J^), and viceversa, i.e. 

(N(u) = 0, 
dE(u)h = VheH <*■ i i _ (14) 

<■ \en 

Proof 

dE\u\-h=—E\u + \h]\ = [ (N'{tu)th,u) + (N{tu),h)dt = 
dX I a=o Jo 

= f (N'(tu)u,th) + (N(tu),h)dt= f ^-(N(tu),th)dt=(N(u),h). 
Jo Jo dt 

So dE[u] = if and only if N(u) = as claimed. 

The finite parameters reduction of the preceeding section can be applied to 
the functional E(u) to obtain a finite parameters variational principle, which is 
equivalent to the infinite dimensional one, i.e. every critical point of E is also a 
critical point for 

E : R m — > R, 

(15) 

This equivalence between the finite parameters variational principle -j-E(fi) = 
and the original Dirichlet problem JT]) is resumed and proved in the following 

Theorem 2 There is a one-to-one correspondence between the critical points jl 
of E(fi) and the solutions u of the Dirichlet problem (l\). 

Proof Plainly, 

dE(fj.) = dE[u] ■ — (g(fi + Tj{n)))dix = 
dfi 

= (N(u)\ g (dfi + fj' \fi)dfi)\ = 

= (¥ m N(u) + Q m N(u) I , g(dfjL) + g(v(n)dA = 

\ Ju=g(p+J7(p)) ^_^^ , ,j 

=0 by JToJi ePm-H eQ m H 

= {¥ m N(g(f, + fi(fM))),g(dfx)). 

We used the facts that (i) the diagonal isomorphism g maps ¥ m H and Q m H into 
themselves and that (ii) V m H and QmH are orthogonal also with respect to the 
L 2 scalar product: 

{ui,Uj) = —{-Lv,i,Uj) = -7- ({ui,Uj} = 0, i* j. 

We can conclude that 

dE(n) = <=> P m N(g(n + ii(n)))=0. 
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4 Quasi- quadratic functions and their properties 

4.1 The deformation Lemma 

Standard cohomology on a differentiable manifold M, 

k i h/r \ {closed k - forms on M} 
{exact fc - forms on M } ' 

is generalized to the relative cohomology by 

H k (M, N) := {q 6 H k (M )| i* a = o} , 

where N is a submanifold (possibly with boundary) of M and t, : N ■-> M is the 
inclusion map. 

It is well known that the compactness hypothesis for the classical calculus of 
variations can be weakened by adopting the Palais-Smale condition on the pair 

(Mjy. 

Definition 1 (Palais— Smale) We will say / : M ->■ R is Palais-Smale, if every 
sequence {a;^} such that |/(£i)l is bounded and /'(xj) -»■ 0, admits a converging 
subsequence. 

Under this condition the critical levels are still compact, even if it could not be 
the case for M. It is possible to associate a critical value 7(a, /) to every non zero 
cohomology class a, i.e., 

Vae#*(M)\{0} 3x e M : f(x) =7(0,/), df(x) = 0. 
Let / := \x e M f(x) < X> denote the sublevel set corresponding to A e R. 
Theorem 3 (minmax) Let / : M ->■ R &e Palais-Smale. Let a e H*{M), a + 0, 

7(a,/):=inf{A|txa*0} (16) 



and i\ '■ f "-* M . Then 



is a critical value for f. 

Proof If c is not a critical value, by Palais-Smale condition, there exists e > such 
that /~ ([c-e, C+e]) does not contain critical points. It is possible to deform f c+e 
in / c_e , by means of the flow of the vector field 

|V/(a;)| 

which could be smoothly extended to all M to the zero vector field outside of a 
neighborhood of / c+e \ / c_e . One can easily check that if tpt is the flow of X(x), 
then 

fWt{x)) = f(x)-t. 

By means of the diffeomorphism ip 2e ■ / c+e -> / c_e , H*(f c+E ) = H*{f c ~ e ), there- 
fore if a ^ 0, it cannot be a = 0. □ 

| /c+e | /c - e 

Note that different cohomology classes could correspond to the same critical value. 
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4.2 Quasi— quadratic functions 

The theory of the Generating Functions Quadratic at Infinity (GFQI) has been 
developed by M. Chaperon, J.C. Sikorav and C. Viterbo, in the context of Hamil- 
tonian dynamics (see 7,8, 15, 20, 21 , 22, 24; ). In the symplectic topology environ- 
ment, D. Theret [22], and later C. Gole [15] . revised carefully results claimed by 
C. Viterbo [25] adopting apparently more general definitions for GFQI. In the 
following we will stick to the Viterbo-Theret format. 

Definition 2 (GFQI) We say that a function S : R ra -» K is quadratic at infinity 
(GFQI) if there exists a non-degenerate quadratic form (Qu,u) and a constant 
K > such that 

\\S(u) - (Qu,u)\\ cl = \\S(u) - (Qu,u)\\ + \\S\u)-2Qu\\ <K (17) 

For this class of functions the Palais-Smale property does hold, and for sufficiently 
large c > 0, the sub-level sets S c := {x : S(x) < c] of a GFQI function S are dif- 
feomorphic to the sub-level sets of the related quadratic form Q, leading as a 
consequence to the equivalence of the relative cohomology classes: 

H*(F- C ,F C ) = H*(Q~ C ,Q C ). 

In this Section we will give an explicit proof of the equivalence S ±c = Q ±c along a 
rather different line of thought than [22] . We start by noticing that all the critical 
points of S are in a compact neighborhood of the origin. Indeed, 



S'(x) = 0, 



\2Qx\<K, m. xzB[0,— — | — ) 

\ 2mm|bpec<y J 



Moreover, GFQI are Palais-Smale. Indeed, let us consider a P-S sequence {^OieN' 
i.e. S(xi) bounded, S'(Xi) -> 0. For i sufficiently large we have ^'(xj)! < K, 

and by quasi-quadraticity, \2Qxi\ < 2K, then Xi e B (0, min ii^ ec0 | )■ Being the 
sequence definitively confined in a compact set, it is possible to extract a converging 
subsequence. We are ready now for proving the main result of this section. 

Theorem 4 For c> sufficiently large, 

i^ diffeo i /. 

S ±c = Q ±c . 

Proof Assume that c is large enough to surpass every critical value of S, and that at 
least one of the eigenvalues of Q is positive, otherwise we easily get S c = Q c = R™. 
Because of (|17p . we have 









• c S c c Q c+K 


c S C+2K 


^Q C 


:+3K 


and also 


















S c 


s S c+K = 


„ qC+M 




and 


Q c = Q c+t 


Therefore 


we 


only need to prove 


















s c ^ Q c+M 







10 
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v(x) = 1 



ip(x) = ■ 



q>(x) = o 







Fig. 2 A continuous Urysohn function ip, s.t. ip = on S c and tp = 1 on Q c 



for some c, M > 0. To this end we define a smooth function / such that f c = S c and 
f c+ = Q c+ , without critical values in [c, c + M] . This can be done by means of 
an Urysohn function ip, 



tp(x) = 



if S(x) < c, 



and setting 



1 if c + M <Q(x), 



(18) 



By doing so the derivative of /: 

|/'| = \[(l-<p)S'+ipQ']-[<p'(Q-S)]\ > \(l - <p)S' + ^'1-^1 -\Q-S\. 



>\Q'\-K 



<K 



\Q'\-K 



is non vanishing if lip is bounded, i.e., \tp < ■ — 4r — . A continuous tp can be defined 
as follows: 



ip(x) 




Q(x)- 



if Q{x) < c, 

if c< Q(x) <c + M, 

\ic+M<Q(x). 



(x)\ \Q'(x)\-K 



(see Figure [2]), indeed |<£>'(:c)| = ■ — j^— - < ■ -^ — wherever ip is differentiable. An 

everywhere smooth Urysohn function can be obtained from tp by mollification. In 
the Appendix we provide an explicit construction and verify that the constraint 
on the derivative still holds. It remains to prove that \Q'\ restricted to {Q = c} is 
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> K, provided that c is sufficiently large. We set v '■= Q'(x) = 2Qx and the problem 
is to minimize ^- for v J -^j-v = c. Lagrange multipliers format gives: 



2 

i 

2 



j-(v,a) = — A I V — — v-c\, 



dT A T 

0= — = v--Q v. 

ov 2 

Therefore, v = 2Qx must be searched among the eigenvectors of Q~ J , or equiva- 
lently x among the eigenvectors m of Q. If x = qu, on the level c, i.e., a = \/vrj 
then we have 

Q'(x) = 2Qx = 2J-^-Qui = 2v^yA~ Ui => |Q'(i)| = 2^7^, 

and the minimum is realized for the smallest positive eigenvalue of Q, showing 
that |Q'(x)| grows at least as y/c on Q = c. □ 



4.3 Existence of a critical point for S 

Clearly, x = is the unique critical point of Q with Morse index equal to q = indQ = 
# {negative eigenvalues of Q). Since 

h , „r „-r- R, if k = a. . 

H k (Q c ,Q c ) = \' (19) 

10, n k + q, 

(see, e.g., [TJ], page 188), we have, for c sufficiently large, 

H*(S c ,S- c ) = H*(Q c ,Q- c )±0 

where also we have used the fact that (relative) cohomology is invariant under 
diffeomorphisms. Thus, there exists a critical value for S between -c and c, corre- 
sponding -by Theorem [3]- to the generator of H* (Q c , Q~ c ). 



5 Quasi-quadratic reduced energy functional 

In this section we will investigate around suitable conditions on F in ([TJ) , in order 
that the associated reduced variational principle being a generating function quasi 
quadratic at infinity, thus admitting at least a critical point. 

We write more precisely (|13p pinpointing the quadratic and the nonquadratic 
part in u: 

E{u) = f (N(tu),u)dt = 
Jo 

= -([ (Uu,u)dt+X f (tu,u)dt\- f (sF(tu),u)dt = 

=:E L+X {u) + E F {u). 



12 
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The reduced functional (|15|) splits consequently, 

E{ji) = E{u{m,v)) = E L+ \u{^)) + E F {u{^)) = E L+ \n) + E F (^), 
and, furthermore, a quadratic form in the /Xj's can be isolated: 



oo I 

j = 1 A J \j = l \ 



(20) 



E L+ \u) = ^ +A ( ff «) = ^ f vj - A | £ ^ I = 



j=m+l A J \j = 

Ai — A _2 



7 = 1 A J j=m+l A J \j = l A j j=m+l A j 



A, -A 



\ 2 



= E 



/*i+ E 






»7j (m) = 



= J B 9 L+A ( / i)+^, +A ( M ). (21) 



Thus we can set: 



Aj-A „2 



E q (n) ■■= T,jLi -a5-Mj (= : (Qm,m>), 

j 

In the general case, there is no reason for i? n q (/i) and -E (/x) to be quadratic. 
With an additional hypothesis on F, it is possible to prove the following theorem. 

Theorem 5 Let the nonlinear functional F : H -> H be a Nemitski operator, 
E(u) = / o u, associated to f : R — ► K. // /, its derivative f and one of its 
primitives f are bounded, then, the reduced functional E(fi) is a GFQI, i.e. 



lE(n)-(Qn,n)\ 



c 1 



< const 



(22) 



where 



1 A-Aj 

A-A, 


O 


O 


A-A m 



\ 

As a result, problem (QJ) admits at least one solution. 

Remark 1 The request on / to be bounded along with its derivative and primitive 
could seem at first overly restrictive to include sufficiently general nonlinearities. 
However, it should be noticed that these restrictions apply to the strictly nonlinear 
part of the operator, and we point out that our class of nonlinearities include Co° 
functions, Schwartz functions, as well zero-mean periodic functions, e.g., sin(x). 



(23) 



Proof Because of (|20[) and (|2ip . we have that 

E(fi) - (Q M , fx) = E{n) - £ 9 L+A (/x) = E^ x (fi) + E F (fi), 

therefore we have to verify that E nq (/1) and E (/1) are bounded along with their 
derivatives. Let K > be the constant such that 



\f\,\f\,\f'\<K. 
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1. \E F \ < cost. Whenever u(x) + 0, 

/ f(tu(x))dt=—- / f(r)dr, 

Jo u(x) Jo 

thus, denoting by Q' = Ix e fi\ u{x) * 0[, 



\e f \ = 



( f eF(tu)dt,u) = f e\—— f" fMdr-uidx 

\Jo I J Q' \ u(x) Jo J 

f ( f^^ f{T)dr\dx -el f f(u(x))d 
Jsr \Jo ) \Jn' 



< eK meas(J?). 



2. |-B^ 9 +A | < cost. We have that 



j=m+l A j 

Being A fixed, while m can be chosen arbitrarily large, it is non restrictive to 
suppose A « Xj, M j > 77i. 



-I oo 



A 



m+l j=m+l 



A, -A 



n 3 

<1 



\Vj(m,fi)\ 2 < IHI 2 

Am+l 



On the other hand, being r\ obtained by means of the fixed point technique, it 
satisfies 

V=(I-\gy 1 eQ m Fg( f i + fj). 
We prove that |(I - A3)" 1 1 < A A "' + ^ A . Indeed, 



-HP 



(l-Xg)x=Ul-—x 3 
(I -A3) ' 



A 

V 

■' ' ) j=m+l 



*i 



=m+l *■ J > J=m+1 



and if we suppose A « A m < Xj , 

Ara+l Aj 



Thus, 



A m +i - A Xj -A j^ 



, , Ajn+1 II ,-,/ / , ~\\M , A m+1 T s 

\n\ ^~ 7£\\F(g(v + v))\\ < reK, 



Am+l — A 



Ara+l — A 



as claimed. 



11 
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To give an estimate on the first derivatives of E and E n g , we start by considering 
fj'(fi) and u'(fj,). Because of the fixed point equation defining fj, we have that 






U(.-w(*».H| ( |)) 



h= 1. 



dfj 



dnh 



\*h A m+ i dfi h ) 
r A m+ i / 1 

T^avaT 



<eK- 



thus, 



dfj 



A, 
eK 



dfj 



dfJ-h 



11'/ 



dfih 
eK A m+ i 
Ai A^ 



L eK \ < eK X m+1 

V Am+i — A/ Ai Am+i — 



1 



An 



A 1-eK- 



eK 



Ai A m+ i - A - eK 



A m+1 -A 

where the rightmost fraction is well defined and positive for m sufficiently large. 
We obtain also an estimate for u , 

u(n) = g{n + f){n)) => u = ff(// + fj'), 

i f\\ ii / ' ~/\ ii , ii ii /-i n~/ii\ ii /ii , ^ /-i n~'ii\ 
\\ u I = ff(M +7 7 )|| ^ Iff I ( 1+ \\V \\),\\u\\<—{l + \\v\\)■ 
3. |-#-.E F | < cost. Recalling that 
I ^m I 

( f F(£u)dt,u) = f f(u(x))dx, 
it is easy to check that 

= \a~ / /(«W)*l! = / /(u(x)) — (x)da 
I o/j, Jn \Jn o/j, 



Oj-L 



4. U^-_E„g < cost. By definition we have 



dE^ x a 



<9/i 9/i 






- / I/I u ' da: < u' if meas(J?). 
-£Vi\ ^ ^ j I A, j' 



then 



9 pL+A, N 



<2 



A 7 - A c?77 7 

A, du . 

J ™ V m+l 

1 

<2- 



A 



■ sup 

m+ 1 j>m 



A, -A 



A, 



F7 -II 7 ?! ^ 



*7 ' ll»7ll 



We can conclude that \\E(fj,) - (Q/j,, /j,)\\ 1 is bounded by a function depending 
only on e, m, K, meas(J2), Ai and A m +i. n 
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6 Optimal mechanical interpretation of the reduction of a continuous 

system 

6.1 A model for the nonlinear elastic string 

In this section we consider a specific continuous problem, whose linear part could 
be obtained as the thermodynamic limit of a sequence of discrete mechanical 
systems. After having applied AZ reduction, by means of a suitable change of 
variables, we restore a discrete mechanical system 

(i) variationally equivalent to the continuum and 

(ii) physically comparable to an element of the sequence involved in the above 
thermodynamic limit. 

Let us consider the one dimensional stationary problem of the nonlinear elastic 
string: 

\au + Tu" = F(u), u:[0,L] 

u(Q)=u(L) = 0. ' 

where u is the vertical displacement of an infinitesimal element of the string, F is, 
as before, a globally Lipschitz nonlinear Nemitski operator, a > or < represents 
a vertical attractive or repelling elastic force, while r > controls the elastic 
attraction among neighboring infinitesimal pieces of the string. The eigensystem 

of ~ T ^ 1 is given by 

(7rfc\ 2 
~r > 1 - fc > 
rr (25) 

uiAx) = \ I sin I — x . 

The above AZ construction leads to the finite variational principle dE(fi) = 0, 
where 

* Xi+o- 2 . + , |F| , (26) 
fj,j, with 



E (A t ) = L^ %> Wlth T^<1 



This exact finite spectral formulation does correspond to the description of a large 
number of more or less realistic finite dynamical systems. The direct reconnaissance 
of (|26p shows a system of N harmonic oscillators, whose positions are described by 
the spectral variables jij, with elastic constants { ; the term E n zn {iJ.) represents 
a further nonlinear coupling, incorporating all the nonlinearities of the system. 

Our aim is to look for a mechanical discrete system, described by genuine phys- 
ical variables, such that, at least when the linear part is concerned, its diagonaliza- 
tion -or spectral representation- is proposing precisely the quadratic part E ITl (/x) 



of (|26|) . Among the large number of possible systems satisfying this requirement, 
we will pick the fitting element from the thermodynamic sequence converging to 
the linear part of the continuous system (|24D . 
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According to this program, we consider the elastic chain of A" + 2 oscillating 
beads (see [131 Chapter 4] and [16]) with fixed extrema: 

Wy i + T(Ay) i = (F(y)) i , i = l,...,N, 
Vo = VN+i = 0, 

where yi represents the i-th vertical displacement, and the finite differences Laplace 
operator: 

i A v)i- - 2 > a= iv7T' l=1 '---< Ar ( 28 ) 

represents the bond among neighboring beads. Each oscillator is constrained to 
move only vertically, is bound to the neighboring beads with a linear spring with 
elastic constant r and also to its rest position at yi = with another spring with 
constant a. In analogy with our continuous problem, we also introduce a nonlinear 
Lipschitz force F depending on the vertical displacements y. 

The role of the variational principle is played by the potential energy, which 
can be written as 

N — N ( \2 i 

u(y) : = z Wi + E V" +l) + uF (y) = \y TA v + uF (y)- ( 29 ) 

Next, we want to exhibit a suitable linear change of variables 

V = y(f-) = ?>, (30) 



such that the potential energy (|29|) assumes the form of (|26|) . More precisely, we 
will choose T in order to make coincide the quadratic parts of (|29|) and (|26p: 



-y(») J Ay(ri = E Hn (v), (31) 

and the remainder representing the overall nonlinear coupling will be defined con- 
sequently: 

U F (y{ti):=Er lin (n). (32) 

Doing so, we would have obtained a discrete system variationally equivalent 
to the nonlinear continuous model and with additionally a strong physical corre- 
spondence, in the sense that the discrete -described by (|29p . and in (|39|) below- 
simulates the microscopic granular structure of the continuum. 

We carry over explicit computations in the following. 



6.2 Equivalence among a chain of springs and the elastic string 

We can compute the normal modes of the linear part of the system by providing 
the full solution set of the homogeneous version of equation (|27p , which is obtained 
via the eigensystem of the finite difference Laplace operator (I28[) : 

Vj = sin I ki 1 , 

- 2 / n \ 4(Af+l) 2 2 ( n \ (33) 

Afc = T— - 1 - cos k =T sin — ; rk . 

a 2 \ N + l ) L 2 \2{N+1) ) 
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It is possible to write the potential energy (J29J) in spectral coordinates £& , obtaining 
an expression analogous to (12(31) : 

N 

k=l 

N 

u(0 ■■= u(y(0) = - £ (<f + A fc ) a 2 + ^ F (C) ( 34 ) 

z fc=i 

= ic T AC + t/ F (C), 

where now J4 is diagonal. Approaching the thermodynamic limit (for JV large), the 
ratio 2IN+1) ~ f° r an y finite k. Therefore we can substitute the sinus with its 
argument in (I33|) : 



- 4(iV+l) 2 . 2 ( JK .\iV-oo 7T 2 a . ,,-x 

At = t— - Sill I — -fe * T^rK = At. 35 

L 2 \2(AT+1) j L 2 V ; 



Therefore, for large JV, the first eigenvalues of the chain approximate the first 
eigenvalues of the continuous elastic string (|25|) . 

If we further apply the following diagonal change of coordinates: 



- + A fe A fc 
2 fc A fc 



(36) 



we obtain a potential energy where the quadratic part now coincides with E (n) 
in ([26)1. 

For the good definition of this transformation "a has to be chosen such that 
the ratio Z. k > for every k. In other words, if a is negative and Xi < -a < Xi+i, 
-a should be placed among the corresponding pair of eigenvalues for the chain: 
Xi < -<r < Ai+i. If a is positive it suffices that also W be positive. 

Finally, composing the two linear coordinate changes (|34ft and (|36|) : 

y = Y( = YZ^i =: Tfi, f i = r- 1 y = Z- 1 Y- 1 y. (37) 



and by substituting in (|26[) , we obtain the finite dimensional non nonlinear energy 
function 

U(y) = E(T- 1 y) (38) 

coinciding with the potential energy of a chain of N + 2 nonlinearly coupled oscil- 
lators with fixed extrema as in (1291): 



N — N { \2 

u(y) = Xn a y> + Lo a — ~2 + u (w)> , TO , 

i=l z i=0 z " (39) 

U F {y) := E nlin {n{y)) = E nlln {T- l y). 
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Fig. 3 Detailed description of the function g. The maximal slope is realized in the central 
interval gz- 



Appendix: A smooth Urysohn function 

We provide an explicit construction of an everywhere differentiable version of the 
Urysohn function in (|18p . 

Lemma 1 Let Q(x) a differentiable function without critical values between c and 
c + M . Then for every e > there always exists a differentiable function tp such 
that 



ifi(x) = 



ifQ{x)<c, 

1 ifc + M <Q(x), 



and furthermore that \(p'(x)\ < ]^j_ for every x. 

Proof By composing Q(x) with a suitable differentiable function g : M. -*■ [0,1], 
strictly increasing in the interval [c, c + M], we will maintain the same sublevel 
sets, and furthermore we are able to control the magnitude of the derivatives. 

In order to define g, we divide the interval [c, c + M] in three subintervals, as 
in Figure [3) 



9(t) = 



9l =Q 

t—c-e 



if t<c, 
g 3 = " M "_2 e ii c + 2e < t < c + M - 2e, 
g 5 = 1 if c + M < t. 



<72 (and similarly (74) should fill the gap with continuous derivatives: 



52(c) =0, 

92(c) =0, 

92(c + 2e)=j^. 



(40) 
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, m=2£ 



'92(c) = 0, 
92(,c+2s) = jj^, 
92(c) = 0, 




/ / m = — 



m-^ 



a = 2e, /? = 



These conditions are met by a unique second order polynomial: 

1 



92(0 : = 

g 4 := c + M ■ 



Ae{M-2e) 
1 



it-cY, 



4e(M-2e) 
We eventually check the bound on the derivatives: 



(t-c-MY 



ip'(x) 



{g[-Q' = if Q(x)<c, 

2e(M-2 £ ) - Q' tfc<Q(x)<c+2e, 



92 -Q' 



93 ■ Q ~ JJZTe 



if c + 2e < Q{x) <c + M-2e, 
94 • Q' = affray Q' Hc+M-2s<Q(x)<c+M, 
[g' 5 -Q' = if c + M<Q{x). 



Therefore we have 

|^(x)| = | 5 '(Q(x))|.|Q'(*)|< 

as wanted. 



1 



(M - 2e) 



IQ'C*)! . 
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